Robustness of Decoherence-Pree States for Charge Qubits under Local Non-uniformity 



O 

o 
o 

Q 



X3 

o 



> 

O 

-I— > 



I 

O 

o 



X 
S3 



Tetsufumi Tanamoto and Shinobu Fujita 
Advanced LSI Technology Laboratory, Toshiba Corporation, Saiwai-ku, Kawasaki 212-8582, Japan 

We analyze the robustness of decoherence-free (DF) subspace and subsystem in charge qubits, 
when difference from the coUective decoherence measurement condition is large (~5%) in the long 
time period, which is applicable for solid-state qubits using as a quantum memory. We solve master 
equations of up to four charge qubits and a detector as a quantum point contact (QPC). We show 
that robustness of DF states is strongly affected by local non-uniformities. We also discuss the 
possible two-qubit logical states by exactly solving the master equations. 



Although decoherence is the largest obstacle for quan- 
tum information processing, a lot of powerful active 
methods for correcting effects of decoherence have been 
discovered^. As to the passive anti-decoherence protec- 
tion, decoherence-free (DF) suhspac^t^ and subsystemm 
are shown to be very useful for collective decoherence en- 
vironment, in which all qubits suffer the same disturbance 
from the environment. Singlet state is the only DF state 
for two qubits and there are two independent DF sub- 
space bases for four qubits, e.g. 
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where |1001)(i234) = |l)i|0)2|0)3|l)4 and so on. The DF 
subsystem starts from three qubits. Experiments have 
succeeded until four qubits in photon system^ and three 
qubits in nuclear magnetic resonance (NMR)i. Bacon 
al^ also showed that, even if there is a symmetry break- 
ing perturbation from the collective environment, which 
is parameterized by a coupling strength 77, the DF sub- 
space is robust in the order of 0{r]) when 77<Cl. 

However, in the case of solid-state qubits, even a single 
qubit is hard to fabricate and the redundancy regarding 
the number of qubits would be a critical issue in con- 
structing a large qubit system. First of all, we could not 
prepare plenty of qubits with mathematically exact size. 
The sizes of Cooper-pair box of Ref.* and GaAs quantum 
dot (QD)^ where coherent oscillation can be observed are 
less than hundreds of nm. The requirement of a few % 
fluctuation between qubits would result in controllability 
of a few nm in fabrication. This would be unrealistic un- 
til future when fabrication process is greatly advanced. 
The fluctuation of sizes would lead to that of interaction 
amplitude between qubits and a measurement appara- 
tus, and that of the applied gate bias, in addition to 
the effect of randomly distributed background trapsiSiii. 
Note that even a roughness of the order of 1 A at inter- 
faces affects current characteristics in advanced LSI tech- 
nologies as shown in Ref.^^ Thus, the non- uniformity of 
solid-state qubits to collective decoherence environment 
is much larger and rather local compared with that of op- 
tical or NMR qubits, and it will be necessary to consider 
the effects of second order 0{rj^) or higher symmetry- 
breaking perturbation. 



(a) four gubits 



(b) three gubits (c) trvo qubits 




FIG. 1: Qubits that use double dot charged states are capac- 
itively coupled to a QPC detector. 



In this paper, we theoretically describe the effect of 
large non-uniformity (~5%) of charge qubit parameters 
on DF states based on time-dependent density matrix 
(DM) equations considering the measurement process by 
detector current. Detection of qubit states induces a 
backaction on the qubit state resulting in a corrupting 
of qubit states. Thus, the measurement is an impor- 
tant interaction with the environment for qubits. The 
charge qubit is a two-level system.^'^ controlled by gate 
electrode si'-S'''^ and constituted from coupled QDs where 
one excess electron is inserted, assuming there is one en- 
ergy level in each QD. The detector discussed here is a 
quantum point contact (QPC) in the tunneling region 
depicted in Fig.l. The position of the excess charge 
affects the QPC current electrically, resulting in detec- 
tion of charged state. Experimentsi^ti^ have success- 
fully proved the high sensitivity of QPC detector current 
/ = GoTvd (Go = 2e^//i ~ 77/xS, T is a transmission 
coefficient Vd is a bias between electrodes) in the tunnel- 
ing region(T < 1). The QPC current induces shot noise 
as the basic and unavoidable noise^^, and the cause of 
decoherence treated here. The purpose of this study is 
to investigate the robustness of many-qubit DF states for 
local large non-uniformity and show the possibility of us- 
ing non-DF states with two-qubit singlet states during 
the QPC measurement. Two-qubit DM is analytically 
solved and perfectly analyzed. 

The Hamiltonian for the combined qubits and the QPC 
is written as H — iJqb-l-i?qpc + -ffint- ^^qb describes the 
interacting iV qubits (Fig.,) : H^b =J2i=i{^i'^ix + eicrtz) 
-+^^^1 Ji.i+iCTizCTi^iz, where ili and are the inter-QD 
tunnel coupling and energy difference (gate bias) within 
each qubit. Here the spin operators are expressed by 
annihilation operators of an electron in the upper and 
lower QDs of each qubit. Ji.i+i is a coupling constant 
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between two nearest qubits, originating from capacitive 
couplings in the QD systemic. | f) and | |) refer to 
the two single-qubit states in which the excess charge 
is localized in the upper and lower dot, respectively. 



s) describes the QPC. Here Cij^s{ci,is) s =Tji is 
the annihilation operator of an electron in the lith (inth) 
level (iLiifi) = l,...,n) of the left(right) electrode. Hint 
is the (capacitive) interaction between the qubits and the 
QPC, that induces dephasing between different eigen- 
states of <Ti2^. Most importantly, it contains the fact 
that localized charge near the QPC increases the energy 
of the system electrostatically, thus affecting the tunnel 
coupling between the left and right electrodes: 
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Note that the case where 5Vi is independent of the qubit 
corresponds to collective environment, thus the DP states 
are realized. Hereafter we neglect the spin dependence of 
V and 5Vi. We assume that the tunneling rate through 
the N qubits, F, is composed of direct series of each N 
tunneling rate near i-ih qubit, P^, such as P^^ = -^7^' 

where Pi is defined as P|^^ = 2iTpLpR{vd/N)\V ± 5V,\'^ 
{pL and pn are the density of states of the electrodes at 
the Fermi surface) depending on the qubit state aiz =H. 
The strength of measurement is parameterized by AF.; 

as P,^^^ = P,o± AP.. We call | ii), | it), | Ti), and 
I IT) as \A) ^ \D) respectively, and four-qubit states 
are written by \AA), \AB) ...\DD). For uniform two 
qubits, Ta = Po(l-C)/2, Fs = Pc = ro(l-C')/2 and 
P£, = P(,( 1+0/2 with C = AP/Po. The DM equations of 
N qubits and detector at zero temperature are derived 
similar to Keii^ as 
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where zi, Z2 = AA, AB, DD for four qubits (256 equa- 
tions) and zi,Z2 — A,B,C,D for two qubits (16 equa- 
tions). JAA^J2i ''i+Jl^J^a, JaB=Y.^ <^t-^4+Jl2-J23, ■■; 

JDD = -J2i ei+Ji2+J23- 9i{zi) and gr{zi) are introduced 
for the sake of notational convenience and determined 
by the relative positions between qubit states. For two 
qubits, 31 (A) = B, 52(A) = C, giiB) = A, g^^B) = D, 
gi{C) = D, g2{C) = A, g^{D) = C, g2{D) = B. Be- 
cause the detector current is / = e^^^^ jjV^pzz, we 
have AT/T ~ AP/Po(= C)^- In the two-qubit case, we 
use entangled Bell basis: \a)= {\A) + \D)) / ^f2, \b)= {\A) 
-\D))IV2, |c)ee m +\C))/V2, \d)^ {\B) -|C))/V2. 

Decoherence rates- From Eq. 0, the dephasing is ex- 
pected to be relevant to the dephasing rate rd(^i, ^2) = 



the rotating 
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^MT^) to eh: 



Tr[p(0)p'(i)] on 

e'S^^''""*p(t)e- 
nate the bonding-antibonding coherent oscillations of free 
qubits (trace is carried out over qubit states). F{t) can 
be expanded in time as F{t) = 1-E,^i(1/"0(Vt^"^)" 
where 1/t(") = {-Tr[p(0)d"/9(0)/dr]}i/" {decoherence 

rates). Using Eq.© for two qubits, l/ri^^ — I/t^^ — 
(l/2)rd(-B,C)~(Po/8)(l-C^)C^772 ^i^gj^ and |C) are 

not symmetric such as the left qubit has a local fluc- 
tuation p(=F)(l-77). Moreover, (1/ri^')^ = (f^i + f^2)^ + 
(ei-e2)' + r|(S,C)/4 and (I/tP)^ = {^,-^2? + ie,- 
€2)^ + rl{B,C)/4:. Thus, the symmetry-breaking terms 
start from 0{rf'^. We have similar expressions for three- 
qubit DF subsystem and four-qubit DP suhspace. For 
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both of which also start from 0{rf') (decoherence rate of 
)(1234) has a similar but more complicated form) 
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The robustness of TV > 3 qubits changes depends on 
which qubit includes local fluctuation. For example, if 



the leftmost qubit of fluctuates 



-^/^l[3] 
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when the rightmost qubit 
fluctuates. In the latter case, there is a symmetry be- 
tween the leftmost and middle qubit and we can say that 
if some symmetry remains, DF states are robust. To 
see the dependence of spatial arrangement of qubits, we 



also consider l^* 
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3 /(1234)= 



7(1423) (hereafter we omit 
the subscript). I^'if') is expected to be more susceptible 
than l^'i'''), because the former is not a product of two 

singlet states as and less symmetric exchanges of 

qubit states are possible when qubit parameters fluctu- 
ate. Thus, the DF states are strongly affected by the 
distribution of local non-uniformities. 

Numerical calculations support these analyses. Here 
we add fluctuations locally to f2.i, and F^ respectively 

to various DF states of A'' = 4 qubits. In case (i), only 
3rd qubit fluctuates as fis ^ ^3(1 — 77), £3 ^ £3(1 — 77) 
and F3 -^F3(l — yy). In case (ii), the 2nd and 3rd qubits 
fluctuate. In case (iii) only 4th qubit fluctuates. Fig. 
12 (a) shows a time-dependent F(t) in a strong mea- 
surement case oi C — 0.6 of 1% fluctuations (?/ = 0.01) 
and Fig 121 (b) shows that of a weak measurement case 
oi C, — 0.2 of 5% fluctuations (7^ = 0.05), both in degen- 
eracy point e = {relaxation decoherence regioiiiS*ii) . 
Although Fig. 121 (a) indicates that the four-qubit DF 
states are fairly robust. Fig. [3 (b) shows that the large 
fluctuation (5%) greatly degrades 1^1'') in case (iii). If 
we take Po=100MHz, its lifetime {F{t) > 0.75) is about 
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500ns, which is much shorter than lowest or- 



/F, 



To see the decoherence effect ex- der estimation of rj Fg 



40/zs. This shows that 



plicitly, we study time-dependent fidelity, F{t) 



higher symmetry-breaking perturbation terms cannot be 
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neglected and there is a case where a product state is 
more preferable than four-qubit DF states. F{t) for |vE'2^') 
is in the same order of j^'^"''), which means that DF states 
composed of many entangled states seem fairly robust 
at the degeneracy point. The three-qubit DF states are 
fairly robust for larger 77, and show susceptibilities to 
measurement strength as the singlet state and other Bell 
states (Figl^a)). In the case of finite bias e {pure dephas- 
ing region) , the effect of which appears from 1/t(2), F{t) 



of \'9^2^) and l^'lj'') seem more susceptible than the singlet 

state and the singlet type l^-J^') (FigEJl. The other Bell 
states and three-qubit DF states are less susceptible than 
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\'^^2^) and I^'Ij''), which would again be due to the higher 
order symmetry-breaking terms. To summarize these re- 
sults, DF states of many qubit {N = 4) are robust for 
most cases but there is a case where even product states 
might be better under large symmetry-breaking fluctua- 
tions [r] ~ 5% ) in the long term. 

Analytical solution for two-qubit case- The probability 
that the unexpected non-uniformities induce symmetry- 
breaking effects becomes higher as the number of qubits 
increases as shown above. Because preparing many solid- 
state qubits is not easy, the redundancy of coding qubits 
is a trade-off against the fabrication difhculty. Thus, 
using the singlet state and one of the non-DF states, 
{paa, Pbb, Pec} in two qubit space would sometimes be a 
realistic solution to construct two logical states \0)l and 
Here we investigate which of {paa, Pbb, Pcc\ is ap- 
propriate for the second basis in two-qubit space. When 
we move to four Bell bases under the conditions that 
the two qubits are identical with no interaction between 
them {Jij =0) and no bias £^ = (collective environment), 
Eqs.Q are divided into the following five groups: 
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Pdd = 

Pbd = Pbd 

Pad = ~2iVLpcd-"f^ Pad 

Pcd = -'^i^Pad 

Pab = -2inpcb+hl^iPba - Pab) 

Pbc ^2iVLpba~l Pbc 

Paa ^2iVL{pac~ Pca)-\l^{Paa 

Pbb ^ [Paa - Pbb) 

Pec ^ -^iflipac - Pea) 

Pac = ^iflipaa ~ Pe ' 



Pbb) 



(5) 
(6) 

(7) 
(8) 

(9) 



-l^Pac 



where 7^ = rd(A,£)) and 7^=rd(A,B). These 
equations can be solved analytically. First, sin- 
glet state Pdd is time-independent (DF state) and 
Pbdit) = pbd(0)e~' ^^*. Pad{t) and pcd{t) behave like 
g-(7^±V(7«)=^-i6n2)t^ ^^^^^^ p^^^^-^ behave like 

g-(7"-+7"±V(7«47")^-i60=)t/2^ thus decay rates of pad{t), 
Pcd{t), Pab{t),and pbe{t) depend on whether 457 > 7^, 
7-^+7^ or 4il < 7^, 7-^-1-7^, respectively. From Eq.®, 
Paait)+pbb{t'yr Pcc{t) is conserved and Pac{t) + Pca{t) = 

(/Oac(O) /9ca(0))e-T«*. ThuS, Paa{t), Pbb{t) , Pec{t) and 
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FIG. 2: Time-dependent fidelity of four-qubit DF states 



I*'*') and ll-^'O) three-qubit DF state of 
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qubit states (singlet, \h) 
der various fluctuations 



= (i-^)r 
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and F, 
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= r 



and ri=^' = (1 



(1- 



and two- 

c), and product state \A)) un- 
(i)ri3 = (l—vi)Q., £3 = riVo and 
n)Q,2 = 0,3 ^ (l-r;)r2, £2 = £3 = n^o 
= (l-r;)r(±\ (iii)!^4 = (1-r;)!^, £4 = ryFo 
~ri)T^^\ (iv)f72 = (1-??)^!, £2 = 77ro and 
r;)r'*'. (iv) is also applied to the two qubits. 



(±) 



(a) 77=0.01 and C=0.6 (strong measurement), (b) r;=0.05 and 
C,—Q.2 (weak measurement). Q, = 2ro, Jij = £i = 0. 



Pac{t) can be analytically obtained by solving three-order 
polynomial equations. In the case of high qubit oscilla- 
tion Vl 3> 7^, 7^ and weak measurement such as C 1, 
we obtain 7^ ^ 47-^ in order of C^. Then we have 
the eigenvalues of time-dependent matrix equations for 
{Paa,Pcc, Pac ~ Pea} in Order of (5 = 7^/(4f]) as -87^, 
— 7^ ± 4rii (1 — 5/3(5^) , which shows that the period of 
qubit oscillation is delayed by the measurement through 
7^ and 7^. Thus we found that the qubit behavior 
strongly depends on the relative magnitude of to 7 
and this is also an important factor in the selection of 
the possible candidate of the logical basis. In Fig.|2b) of 
f7 = 2F > 7^ - 0.2F, F(i)s of \h) and \d) are larger than 
those of |c) (and \a)). In RefJ^, we argued the coherent 
oscillation of the DM without the detector and showed 
that the DM for |5) and \d) are time-independent at e = 0. 
This indicates that \b) is the candidate in the 3> 75, 7_d 
region, reflecting that \b) becomes an unpolarized triplet 
suffering less degradation from the repulsive Coulomb in- 
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FIG. 3: Fidelity of four-qubit DF states at t = 50ro^in the 
case of (iii) as a function of non- uniformity, = 2ro, Jij — 
and C=0.2. 
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FIG. 4: Linear entropy Sl and concurrence C plane under 
the measurement during = t = lOOFo for two qubit states 
\b) and |c) (e — 0). In weak measurement case of = 0.2, C 
does not degrade to 0. All data start closely to Werner state. 



teraction of the detector. On the other hand, in the = 
case which can be in the limit of <C 7^,7^, Eq. Q 
is easily solved even with a finite bias e and pcc is found 
to be a time- independent unpolarized triplet, being the 
candidate of logical state. 

Next, we compare pbb and pcc in the purity plane, that 
is a relation between linear entropy and concurrence^'^. 



Linear entropy Sl = 4/3(1 — Tr(p^)) expresses purity of 
qubits, ranging from (pure state) and 1 (maximally- 
mixed state). From Eq.©, Sl = 4/3[l - pl^ - pg^ - 
pI^ - 2|pa,|2], and we have dSL/dt= {2/3h^{pbb- Paa)^ 
-|-(8/3)7^|pacP- Concurrence C, which is a measure of 
entanglement, is also given from Eq.®. Starting from 
Pbb{0) — 1, at i ^ in the case of ^ 'fn — 47_b, we 
have C - (-l-|-4e-37^t)/3 ( dC/dt^ ~j^pbb/2) and 
Sl 8/9(1-6-67^*). Then dC/dSL —3/4 near t - 0. 
If we start Pcc(O) = 1 in the case of <C 7_d = 47s, we 
have C ~ (-H-4e-(2/3)f^*) and 5^ - 8/9(l-e-(4/3)ot), 

and thus dC/dSL ^^-3/4. If we check the Werner state 
7iii|^) (^1 + (1 — 7'!ti)/4/2 ® I2 where I2 is 2x2 unit matrix 
(1 > 7i„ > 0)22*21, we obtain dC^""^ / dS^^'' = -3/4 at 
7u) = 1. Thus both pbb and pcc coincide with the Werner 
state at t ~ 0, which shows that the two states are good 
entangled states, because the Werner state is a mixture 
of the maximally entangled state2£. Figure 01 shows that 
both Pbb and Pcc evolve close to the Werner state. Thus, 
the two states behave similarly in the purity plane, and 
are equal candidates. 

The noise spectrum S{ijj) of the QPC without qubits 
is given by ^(a;) = b^F/tt (white noise) in the present 
model, thus, the shot noise affects qubit states in full fre- 
quency domain. Astafiev et alm^ experimentally showed 
that the main causes of the noise in the Josephson qubits 
are / noise and the background charge noises or 1/f 
noise, which we do not include. These noises would lo- 
cally affect qubits and degrade the robustness of the DF 
states more than discussed here. 

In conclusion, we have solved master equations of many 
qubits and QPC detector, and discuss the robustness of 
DF states under the large non-uniformities (~ 5%) dur- 
ing the long time period. Two-qubit non-DF states are 
shown to be one solution in constructing logical qubits. 
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